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(57) Abstract 

The invention relates to a method for low-fuel, computer-assisted control of any number of thrusts mounted in any given arrangement 
on a spacecraft, which solves the problem of linear optimization by reducing computational effort in a very short time. The method is based 
on the dual simplex algorithm and suitable provision of a dual admission start table. The dual admission start table is set up on the basis of 
a previously calculated optimum table and a transformation of the actual force/power vector pertaining to the chosen optimum table. Given 
processing steps of the dual simplex algorithm are skipped totally or are substituted by processing steps which require substantially less 
computational effort. Finally, a reduced table as opposed to the simplex table is used. Said reduced table consists solely of an identification 
number for selecting a thruster from the entire number of thrusters and the vector of the control values of the thrusters which are contained 
in said selection. Reduction in computational effort is achieved by carrying out given processing steps in advance and by providing the 
corresponding results in the form of data fields and index or identification number matrixes in the working memory of the computer. Said 
results from the initialization phase of the method need only to be addressed and retrieved in order to calculate thruster control. 



(57) Zusammenfassung 

Es soil ein Verfahren zur insgesamt treibstoffminimalen, rechnergestutzten Ansteuerung einer beliebig festgelegten Anzahl und beliebig 
festgelegten Anordnung von Diisen an einem Raumfahrzeug angegeben werden, das dieses lineare Optimierungsproblem durch Reduzierung 
des Rechenaufwands in mOglichst kurzer Zeit !6st. Das Verfahren basiert auf dem dualen Simplexalgorithmus und der geeigneten 
Bereitstellung eines dual zulassigen Starttableaus. Das dual zulassige Starttableau wird aus einem vorab berechneten Optimaltableau 
und der zu diesem gewahlten Starttableau dazugehorigen Transformation des aktuellen Krafte-Momenten-Vektors aufgebaut. Bestimmte 
Verarbeitungssschritte des dualen Simplexalgorithmus entfallen vollstandig bzw. werden durch Verarbeitungsschritte substituiert, die 
wesentlich weniger Rechenaufwand benStigen. SchlieBlich wird nicht mit dem Simplextableau, sondern mit einem reduzierten Tableau 
gearbeitet. Dieses reduzierte Tableau besteht lediglich aus einer Identifikationsnummer fur eine DUsenauswahl aus der Gesamtzahl der Diisen 
und dem Vektor der Ansteuerwerte fQr diejenigen Diisen, die in dieser DOsenauswahl enthalten sind. Die Reduzierung des Rechenaufwands 
wird dadurch erreicht, daB bestimmte Verarbeitungsschritte vorab durchgefuhrt werden und die entsprechenden Ergebnisse dann in Form 
von Datenfeldem und Indizes- bzw. Identifikationsnummemmatrizen im Arbeitsspeicher des Rechners bereitliegen. Diese Ergebnisse des 
Initial isierungsteils des Verfahrens brauchen dann zur Berechnung der DUsenansteuerung nur noch adressiert und abgerufen zu werden. 



LEDIGUCH ZUR INFORMATION 

Codes zur Identifiziening von PCT-Vertragsstaaten auf den Kopfbflgen der Schriften, die internationale Anmeldungen gemass dem 
PCT veroffentlichen. 



AL 


Albanien 


ES 


Spanien 


LS 


Lesotho 


SI 


Slowenien 


AM 


Armcnicn 


FI 


Finnland 


LT 


Litauen 


SK 


Slowakei 


AT 


Osterreich 


FR 


Frankreich 


LU 


Luxemburg 


SN 


Senegal 


AU 


Australien 


GA 


Gabun 


LV 


Lett land 


sz 


Swasiland 


AZ 


Aserbaidschan 


GB 


Vereinigtes Konigrcich 


MC 


Monaco 


TD 


Tschad 


BA 


Bosnien-Herzegowina 


GE 


Georgien 


MD 


Republik Moldau 


TG 


Togo 


BB 


Barbados 


GH 


Ghana 


MG 


Madagaskar 


TJ 


Tadschikistan 


BE 


Belgien 


GN 


Guinea 


MK 


Die ehemalige jugoslawische 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Griechenland 




Republik Mazedonien 


TR 


TUrkei 


BG 


Bulgaricn 


HU 


Ungam 


ML 


Mali 


TT 


Trinidad und Tobago 


BJ 


Benin 


IE 


Irland 


MN 


Mongolei 


UA 


Ukraine 


BR 


Brasilien 


IL 


Israel 


MR 


Mauretanien 


UG 


Uganda 


BY 


Belarus 


IS 


Island 


MW 


Malawi 


US 


Vereinigte Staaten von 


CA 


Kanada 


IT 


I tali en 


MX 


Mexiko 




Amerika 


CF 


Zentralafrikanische Republik 


JP 


Japan 


NE 


Niger 


UZ 


Usbekistan 


CG 


Kongo 


KE 


Kenia 


NL 


Niederlande 


VN 


Vietnam 


CH 


Schweiz 


KG 


Kirgisistan 


NO 


Norwegen 


YU 


Jugoslawien 


CI 


Cote d'lvoire 


KP 


Demokratische Volksrepublik 


NZ 


Neuseeland 


zw 


Zimbabwe 


CM 


Kamenm 




Korea 


PL 


Polen 






CN 


China 


KR 


Republik Korea 


PT 


Portugal 






CU 


Kuba 


KZ 


Kasachstan 


RO 


Rumanien 






cz 


Tschechische Republik 


LC 


St. Lucia 


RU 


Russische Federation 






DE 


Deutschland 


U 


Liechtenstein 


SD 


Sudan 






DK 


Danemark 


LK 


Sri Lanka 


SE 


Schweden 






EE 


Estland 


LR 


Liberia 


SG 


Singapur 







WO 98/49058 * ' PCT/EP98/01595 



Verfahren zur treibs tof f minimalen , rechnerges tutzten Ansteue- 
rung fur beliebig an einem Raumf ahrzeuq angeordnete Diisen 



Die Erfindung bezieht sich auf ein Verfahren zur insgesamt 
treibs toff minimalen, rechnerges tutzten Ansteuerung einer 
beliebig festgelegten Anzahl und beliebig festgelegten 
Anordnung von Diisen an einem Raumf ahrzeug . Die Diisen konnen 
dabei diskret Oder auch kont inuierlich ansteuerbar sein. 

Fur die Ausiibung von aufleren Kraften und Momenten an einem 
Raumfahrzeug ist der Einsatz von Dusen in vielen Situationen 
notwendig. Da der Einsatz der Diisen Treibstoff verbraucht und 
dieser unmittelbar Einflufl auf das Startgewicht und die 
Lebensdauer von Raumf ahrzeugen und insbesondere von Satelliten 
hat, kommt der Dusenans teuerung mit in der Summe minimalem 
Treibstoff verbrauch eine wichtige Bedeutung zu. Dusenanordnun- 
gen, die prinzipiell einen niedrigen Treibs tof f verbrauch 
ermoglichen, eine hohere oder gleiche Redundanz bei kleinerer 
Dusenanzahl aufweisen und/oder durch die Nutzlast am Raumfahr- 
zeug beeinflufit sind, sind oft unregelmaBige und schief winklige 
Diisenanordnungen. Aus diesen Griinden ist die Forderung nach der 
treibstoff minimalen DUsenansteuerung , insbesondere auch fur 
beliebige ( schief winklige ) feste Diisenanordnungen an Raumfahr- 
zeugen von besonderer Relevanz . 

Es ist bekannt, daG sich das Problem der treibs tof f minimalen 
DUsenansteuerung auf ein lineares Optimierungsproblem zuriick- 
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fiihren laGt, fur das der Simplexalgori thmus ein geeignetes und 
effizientes Losungsverf ahren darstellt (siehe zum Beispiel 
P.J. Wiktor "Minimum Control Authority Plot: A Tool for 
Designing Thruster Systems" In: Journal of Guidance, Control, 
5 and Dynamics, Vol. 17, No. 5, Sept. -Oct. 1994, pp. 4-5). Der 

Einsatz des Simplexalgori thmus scheitert jedoch haufig an der 
bendtigten Rechenzeit. So konnen oft die Echtzeitanf orderungen 
zum Beispiel fiir den Einsatz des Simplexalgorithmus innerhalb 
des Regeltaktes des Lage- und Pos itionsregelkreises eines 
10 Raumfahrzeuges nicht erfiillt werden, weil die besonders in 

hohen Orbits einsetzbaren Rechner den Simplexalgorithmus nicht 
ausreichend schnell abarbeiten konnen. 



Eine Vielzahl der sich heute im Einsatz befindenden Raumfahr- 

15 zeuge (Satelliten) besitzen eine Diisenanordnung, bei der jede 

Diise zum kdrperfesten Koordinatensystem des Raumf ahrzeugs zu 
einer Achse parallel und damit zu den zwei anderen Achsen 
orthogonal angeordnet ist. Diese orthogonale und 
achsenparallele Anordnung hat zur Folge, daB jede Diise nur 

20 Krafte in einem Richtungs sinn einer Achsenrichtung und Momente 

nur in einem Drehsinn urn die beiden anderen Achsen erzeugen , 
kann. Der Vorteil dieser Anordnung ist, dafi die Steuerung eines 
jeden Freiheitsgrades des Raumfahrzeuges als starrem Korper 
durch den Einsatz von jeweils nur zwei Diisen moglich ist, ohne 

25 andere Freiheitsgrade zu beeinf lussen . Die Ansteuerung der 

Diisen ist fur diese Diisenanordnung sehr einfach. Der Nachteil 
dieser Diisenanordnung besteht darin, dafl mindestens zwolf Diisen 
notwendig sind, urn alle Freiheitsgrade gleichzeitig und 
beliebig steuern zu konnen. MuB der Ausfall von zwei beliebigen 

30 Diisen tolerierbar sein, sind 36 Dusen erf orderlich . 

Demgegeniiber liegt der Vorteil einer geeigneten schiefwink- 
ligen, also nicht orthogonalen und achsenparallelen Diisen- 
anordnung im wesentlichen in folgenden drei Punkten: 



^5 



1. Bei frei wahlbarer fester Diisenanordnung sind weniger als 

12 Dusen erf orderlich. 

2. Fiir die Tolerierung des Ausfalls von zwei beliebigen Dusen 
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sind erheblich weniger als 36 Diisen erf orderlich . 
3. Bex treibstof f minimaler Ansteuerung der Dusen ist ein 
geringerer Treibstof fverbrauch moglich. 

5 Der Nachteil dieser Diisenanordnung besteht darin, dan die 

treibstof fminimale Diisenansteuerung die Losung eines linearen 
Optimierungsproblems erfordert und der Rechenauf wand fur die 
Losung dieses Problems mit dem Simplexalgorithmus oft zu hoch 
ist. Aus diesem Grunde sind andere Verfahren entwickelt worden, 
10 die weniger Rechenzeit benotigen und durch die eine mehr oder 

weniger treibstof fminimierte , aber keine insgesamt treibstoff- 
minimale Diisenansteuerung erreicht werden kann. 

Beim sogenannten "Linear Programming" wird das Diisenansteue- 
15 rungsproblem als lineares Optimierungsproblem aufgefaflt und mit 

einschlagigen Verfahren der linearen Optimierung, wie zum 
Beispiel den Simplexalgorithmus, gelost. Der schon erwahnte 
Nachteil, der den Einsatz dieses Verfahrens in vielen Fallen 
unmoglich macht, ist der zu hone Rechenauf wand . 

20 

Beim "Table Look-up" -Verfahren werden alle treibstof fminimalen 
Dusenanordnungen fur Einheitskraf te und -momente bezuglich 
einer Achse des korperfesten Koordinatensys terns des Raumfahr- 
zeugs berechnet ( Losen des linearen Optimierungsproblems) und 

25 in einer Tabelle abgelegt. Sind im Einsatz Krafte und Momente 

bezuglich mehrerer Achsen gleichzeitig auszuiiben, werden die 
entsprechenden Dusenans teuerungen bezuglich Kraft Oder Moment 
und fiir jede geforderte Achse und dem geforderten Richtungs- 
oder Drehsinn aus der Tabelle ausgelesen und mit dem entspre- 

30 chenden Betragswert der Kraft- bzw. Momentenkomponente in 

dieser Richtung multipliziert . Den Ansteuerwert fiir jede Duse 
erhalt man durch die Addition der einzelnen Werte fiir die 
jeweilige Dtise. Dieses Verfahren ist sehr einfach und sehr 
schnell, fuhrt im allgemeinen aber zu nicht treibstof fminimalen 

35 Dusenans teuerungen. 

Bei verfeinerten und erweiterten "Table Look-up" -Verfahren 
handelt es sich urn Verfahren, die eine Look-up-Tabelle anlegen, 
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die Eintrage fur zusatzlich eingefuhrte Achsen enthalten. Durch 
diese Verfahren konnen t reibs tof f sparende Diisenansteuerungen 
ermittelt werden. Zugleich nimmt aber die Obersichtlichkeit des 
Verfahrens ab und wegen der zunehmenden GroBe der Tabellen der 
Rechenauf wand erheblich zu. 

Andere Verfahren beruhen auf der Auswertung und Verbesserung 
des Ergebnisses, das mit dem "Table Look-up" -Verfahren 
ermittelt wird. Das Ziel dieser Auswertung besteht darin, 
besonders t reibs tof f verschwendende Diisenansteuerungen zu 
verhindern bzw. geeignet zu korrigieren. Hier soli exemplarisch 
ein Verfahren, das auf einem Fuzzy-Ansatz beruht, genannt 
werden (vgl. T.Suzuki, K.Yasuda, S.Yoshikawa, K.Yamada, 
N.Yoshida "An Application of Fuzzy Algorithm to Thruster 
Control Systems of Spacecraft" In: ISTS 94-C-16, 19th Inter- 
national Symposium on Space Technology and Science, Yokohama, 
Japan, May 15-24, 1994, pp. 1-8). 

Eine Treibstof f minimierung der Dusenans teuerung insgesamt hat 
zum Ziel eine solche Ansteuerung der Diisen zu ermitteln, die 

1. die geforderten Krafte und Momente am Raumfahrzeug ausiibt, 

2. fur jede Diise einen nichtnegativen Ansteuerwert liefert, 
und 

3. die Summe aller Ans teuerwerte der Diisen minimiert. 

Dies stellt - wie erwahnt - ein lineares Optimierungsproblem 
dar. Da das Losen des Problems in der Regel zu zeitauf wendig 
ist, kann die Losung des quadra tis chen Opt imierungsproblems als 
Ersatzproblem durchgefiihrt werden. 

Die Ldsung des Ersatzproblems fiihrt zur Ermittlung derjenigen 
Dusenans teuerung, die 

1. die geforderten Krafte und Momente am Raumfahrzeug ausiibt, 

2. fiir jede Diise einen nichtnegativen Ansteuerwert liefert, 

und 
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3. die Summe aller quadratischen Ans teuerwerte der Diisen 
minimiert . 

Dieses quadratische Optimierungsproblem kann uber die Bereit- 
stellung der Moore-Penrose-invertierten Sys temmatrix der 
Dusenanordnung am Raumf ahrzeug und eines Of f setvektors , der mit 
Hilfe der singularen Wertezerlegung berechnet wird (vgl. z.B. 
B.Noble, I,W. Daniel "Applied linear algebra" Prentice-Hall 
International, INC • (1988), pp. 338 - 350), geldst werden. Die 
Losung des quadratischen Opt imierungsproblems kommt dann als 
Produkt der Moor-Penrose-invertierten Systemmatrix mit dem 
kommandierten Kraf te-Momenten-Vektor und der Addition des 
Of f setvektors zustande. Die Addition des Of f setvektors bewirkt, 
dafi im Losungsvektor nur nicht negative Ans teuerwerte 
auftreten. Die Berechnung des Of f setvektors wird jedoch dann 
numerisch immer aufwendiger, wenn die Zahl der anzus teuernden 
Diisen um zwei Oder mehr uber der Dimension des kommandierten 
Kraf te-Momenten-Vektors liegt. Im allgemeinen ist die Losung 
nicht treibstoff -minimal , fallt aber mit der Losung des 
linearen Optimierungs-problems zusammen, wenn die Zahl der 
anzusteuernden Diisen nur um eins grofter ist als die Dimension 
des kommandierten Kraf te-Momenten-Vektors (vergleiche. z.B. 

DE 195 19 732 A1 ) . Weiterhin ist nachteilig, dafl langere 
Einschaltzeiten ttberproportional "bestraft" werden und daher 
seltener in den Ans teuerungen vorkommen. Fur den Betrieb der 
Diisen und ihrem Wirkungsgrad ist es jedoch giinstiger, sie mit 
langeren Einschaltzeiten zu betreiben. 

Der Nachteil bestehender Verfahren besteht mithin darin, daft 
sie entweder nicht das lineare Optimierungsproblem losen und 
damit im allgemeinen nicht die treibs tof f minimale Ansteuerung 
der Diisen berei ts tellen oder aber, wenn bei den Ansteuerver- 
fahren das lineare Optimierungsproblem geldst wird und der 
Simplexalgorithmus ausgefiihrt wird, dieses ublicherwei se zu 
rechenauf wendig ist. 



Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der 
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eingangs angegebenen Art bereitzustellen, mit dem das lineare 
Optimierungsproblem fiir die treibstof f minimale Ansteuerung 
beliebig fester Dusenanordnungen zur gleichzei tigen Reali- 
sierung von Kr&ften und Momenten an einem Raumf ahrzeug in 
signifikant kiirzerer Zeit gelost werden kann, als es bisher 
mdglich ist. 

Diese Aufgabe wird gemafi der Erfindung durch die im Anspruch 1 
gekennzeichneten Merkmale gelost. 



Das Verfahren nach der Erfindung basiert also vorteilhaft auf 
einem, dem Simplexalgorithmus angelehnten Verfahren, das jedoch 
aufgrund der konsequenten Ausnutzung der Besonderheiten bei der 
Dusenansteuerung bei Raumf ahrzeugen und der Aufgabe der Verwen- 
15 dung des iiblichen Simplextableaus zu einer Rechenzeitverkurzung 

von ca. einer Groflenordnung fuhrt. 

Vorteilhafte Ausgestaltungen und Weiterf uhrungen des Verfahrens 
nach der Erfindung sind in den dem Anspruch 1 nachf olgenden Pa- 
20 tentanspriichen gekennzeichnet . 

Das Verfahren nach der Erfindung basiert auf folgenden Grund- 
lagen : 

a) Das bei der Dusenansteuerung zu losende lineare Optimie- 

25 rungsproblem kann wie folgt auf geschrieben werden: 

Gesucht ist der Ansteuervektor k ^, der die Ans teuerwerte fiir 

—opt' 

alle Dusen enthalt und folgende Bedingungen erfullen mufl: 

1. Durch k . miissen die geforderten Krafte und Momente b, 

—opt — kom 

30 am Raumf ahrzeug ausgeiibt werden. Das heifit, die Glei- 

chung A . k , = b. muli erfiillt werden, wobei A die 
— sys —opt — kom ' — sys 

Systemmatrix ist, in der die drei Kraft- und die drei 
Momentenkomponenten fiir jede einzelne Duse in einer Spalte 
dieser Matrix eingetragen sind. Die Systemmatrix hat 6 Zeilen 
35 wegen der 6 moglichen Freiheitsgrade in der Bewegung des 

Raumf ahrzeugs , namlich 3 Krafte in drei zueinander orthogonalen 
Achsen und entsprechend 3 Momente urn diese Achsen und bei 
D Dusen D Spalten. 
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2 - iS:^^. = Ck . k , . . . k _ ] darf nur nichtnega- 

opt opt 1 opt j opt D * 

tive Werte enthalten, 

d * h - : k opt j * °' ^ = 1 (1 )D * 
3. die Summe aller Ans teuerwerte muft minimal sein 

( Zielf unktion) , 



d.h. 



D 

E k ^ .= Min ! 
j-1 ° Pt J 
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15 



b) Ein Simplextableau besteht aus einem Spaltenvektor der 

Basisvariablenindizes \ einem Zeilenvektor der Nichtbasis- 

variablenindizes % der Koeff izientenmatrix A^ dem 

— nb' —tab' 

Zeilenvektor der auf die Nichtbasis variablen umgerechneten 

Zielf unktionskoef f izi^nten £ tab und dem Spaltenvektor der 

Basisvariablenwerte b , : 

—tab 



20 



25 







. . • 3 

A nbji+j 








a U 


• • • a xj 


"' a lm 


b X 




a kl 


a kj 


• 
• 






a nl 


• 










... p 







Im Verfahren nach der Erfindung kommt das Simplextableau nicht 
30 in dieser vollen Form zum Einsatz. 

c) Fur die Umsetzung eines beliebigen 6-dimensionalen Krafte- 
Momenten-Vektors sind maximal 6 Diisen einer vorhandenen festen 
Diisenanordnung erforderlich Oder der kommandierte Krafte- 
Momenten-Vektor ist mit der vorhandenen Diisenanordnung nicht 

>5 realisierbar . 

d) Fur jede. Diisenauswahl I von 6 aus D Diisen existiert ein 
Basisvariablenindizesvektor X (I) . Die Nummern oder Indizes 

— D 

(1,2,...,D) der in der Auswahl I enthaltenen Diisen werden im 
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Basisvariablenindizesvektor X b (D zusammengef aBt und die 

Ansteuerwerte dieser Diisen sind die Werte der Basi s variablen 

des Simplextableaus . Die Werte der Basisvariablen sind 

5 entsprechend der Reihenfolge der Dusenindizes in% b (I) im 

BasisvariablenwerteveJctor b^ .(X^, b ) abgelegt. Zu iedem 

—tab — b'— com 3 3 J 

Basisvariablenindizesvektor X^CI) gehdrt eine Basis- 
matrix Die Basismatrix wird gebildet, indem man 

diejenigen Spalten der Systemmatrix A , die zu den 

* —sys' 

10 ausgewahl ten Diisen gehdren, zu einer neuen Matrix spaltenweise 

zusammenstellt . Die Reihenfolge der Spalten in der Basismatrix 
entspricht der Reihenfolge der Indizes der Diisen im 
dazugehorigen Basisvariablenindizesvektor X (I) . Die 
Nichtbasismatrix B nt) (I) wird aus denjenigen Spalten der 

15 Systemmatrix zusammenges tellt , die zu den nicht in der Auswahl 

I enthaltenen Diisen gehdren. Ihre Indizes sind im Nichtbasis- 
variablenindizesvektor X^* 1 ) zusammengef afit . Die Reihenfolge 
der Spalten in der Nichtbasismatrix entspricht der Reihenfolge 
der Indizes der Diisen im dazugehorigen Nichtbasisvariablen- 

20 indizesvektor. Die Ansteuerwerte dieser Diisen sind Null. 

e) 1st fur einen kommandierten Kraf te-Momenten-Vektor b, 

— Kom 

das Optimaltableau X w l , X . ^ , a^ ^ p x u und 

—b.opt ' —nb.opt' —tab. opt' -^tab.opt 

—tab opt 9 efunden ' dann sind die Ansteuerwerte der Diisen, die 

durch die Indizes in "X adressiert werden, gleich Null und 

—nb.opt 

25 diejenigen Diisen, die den Basisvariablen entsprechen, also 

durch die Indizes in X u adressiert werden, enthalten als 

—b.opt 

Ansteuerwert den Wert aus der entsprechenden Zeile des opti- 

malen Basisvariablenwertevektors b^ . t . Der Basisvariablen- 

— tan . opt 

indizesvektor X. gibt die optimale Diisenauswahl an. 

— o.opt 

30 f) Fiir den Fall der gleichzeitigen Ausiibung von alien drei 

Kraf tkomponenten und alien drei Momentenkomponenten ist der 
Einsatz von maximal 6 der vorhandenen Diisen erf orderlich. Es 
gibt genau 



D! 



A(C 6 ) 



35 ~ 6! (D-6)! 

Moglichkeiten ( Kombinationen ohne Wiederholung) , 6 aus D Diisen 
auszuwahlen. Werden die Indizes in den beiden Vektoren der 
Basisvariablenindizes und Nichtbasis var iablenindizes , "X. und 
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— nb' stets aufsteigend sortiert, so gibt es jeweils ge- 

nau A(c£) dieser Vektoren. Fiir eine feste Dusenanordnung gibt 

es (nur) eine bestimmte Anzahl von Kombinat ionen von Diisen, die 

die Diisenauswahl fur eine Losung des linearen Opt imierungs- 

5 problems sein kdnnen. Diese Anzahl ist im allgemeinen deutlich 

kleiner als die Anzahl A(C^). 

6 

Eine wesentliche Grundlage des Verfahrens besteht darin, dafi 

stets mit sortierten Basis- und Nichtbasisvariablenindizes- 

10 vektoren gearbeitet wird. Im weiteren sind die Indizes oder 

Nummern der Diisen in den Basis- und Nichtbasis variablenindizes- 

vektoren z.B. stets der Grofie nach sortiert. Auf diese Weise 

wird es moglich, dali einer Diisenauswahl genau ein Indizes- 

vektor, genau eine Koef f izientenmat rix A = [a .], genau 

t a d tao k.j 

15 ein Koef f izientenvektpr der auf die Nichtbas i svariablen 

umgerechneten Zielfunktion 2 tab und genau eine Basismatrix B fa 
und Nichtbasismatrix B^ zugeordnet werden. 

Eine weitere wesentliche Grundlage des Verfahrens besteht 

20 darin, daG mit einem reduzierten Tableau gearbeitet wird. Das 

reduzierte Tableau besteht lediglich aus der Diis enauswahl- 

identif ikationsnummer I und dem Basisvariablenwertevektor b . . • 

~taD 

Fiir das Verfahren und die Arbeit mit dem reduzierten Tableau 
ist neben der Bereits tellung des Feldes der Koef f izienten- 
25 matrizen aller Optimaltableaus und der Pivot spal tenindizes- 

matrix auch die Berei ts tellung der Dii senauswahl identif ika- 
tionsnummernmatrix, der Zeilenaus tauschindizesmatrix und der 
Basisindizesmatrix erf orderlich . 

30 Das Verfahren besteht aus zwei Bes tandteilen . Der erste, der 

Initialisierungsteil , dient zum einen der Be reit stellung des 
Feldes der Koef f izientenmatrizen aller Optimaltableaus 

£ = 4] = EM' wobei [*** 4 = ^ (/) = ^ (/) ist - 



35 



Die Dusenauswahlident if ikationsnummer lauft von 1 bis 

M ( I = 1(1 )M ), wobei M die Anzahl aller Dusenauswahlmoglich- 

keiten ist, die fur einen beliebigen Kraf te-Momenten-Vektor 
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— kom Und die vorlie 9 ende beliebige und unveranderliche 
Dusenanordnung optimal sind. Eine Koef f izientenmatrix eines 
Optimaltableaus liegt im Sinne dieses Verfahrens dann vor, wenn 
alle Elemente ira zur Koef f izientenmatrix gehorenden Vektor der 
5 auf die Nichtbas isvariablen umgerechneten 

Zielf unkt ions koef fizient en 

£ tab = ^tab'- tab* nicht negativ sind. Fur die Elemente der 
Koef f izientenmatrix mufi daher gelten: 

6 

10 1 - £ a.. . = p . * 0 fur j = 1(1)6. 

k==1 tab k j r tab j J 

Zum anderen werden im Initialisierungsteil des Verfahrens 
f olgende Adreflmat rizen : die Pivot spal ten indizesma trix 
P - tp Ik ] # die Dusenauswahlidentif ikat ionsnummernma- 

15 trix U - C u Ik 3, die Zeilenaus tauschindizesma tr ix T = ^j}^ und 

die Basisindizesmatrix L = [1 ] bereitges tell t , damit auf die 
im zweiten Teil des Verfahrens benotigten Daten besonders 
schnell zugegriffen werden kann und bei der Durchfuhrung der 
Basisaustauschschritte die Arbeit mit dem reduzierten Tableau 

20 moglich wird. Durch die Dusenauswahlidentif ikationsnum- 

mer I = 1,2, M wird jeder Diisenauswahl , die optimal sei-n 

kann, und damit auch jeder Koeff izientenmatrix eines Optimal- 
tableaus eine Identif ikationsnummer zugewiesen. 

25 In der Pivot spaltenindizesmat r ix P = Cp Ik 3 wird fur jede Zeile 

von jeder in F abgelegten Koef f izientenmatrix der Index der 

dazugehdrigen Pivotspalte s € [0,1,2,3, , (D - 6)} abgelegt. 

Fur eine Koef f iz ientenmatrix A tab und eine gewahlte Pivotzei- 
le z ist die Pivotspalte s diejenige, fur die gilt: 



Minimum Pj 



a 



num. y: D 

■< <0 pj = i5' mi,/ - l(l)(z> - 6 »- 



D,h. die zu einer Zeile einer Koef f izientenmatrix gehorige 
Pivotspalte ist diejenige Spalte, in der das Mat rixelement der 
35 Koef f izientenmatrix in dieser Zeile negativ ist und der 

Quotient aus dem Element des Vektors der auf die 

Nichtbasisvariablen umgerechneten Zielf unktionskoef f izienten in 
der Spalte und dem Betrag des ent sprechenden Mat rixelements aus 
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der Koef f iz ientenmatrix am kleinsten ist. 

Mit Hilfe der Pivotspaltenindiz esmatrix ist fur eine Diisen- 

auswahl I, fur die in F eindeutig die dazugehorige 

Koef fizientenmatrix A tab (I) abgelegt ist, und gewahlter 

Pivotzeile z die Angabe der ent sprechenden Pivotspalte mit 

s(A tab (I),z) = P Iz au s der Pivot spaltenindizesmatrix moglich, 

ohne daG dazu die Koef fizientenmatrix A, . ( I ) Oder der Vektor 

—tab 

der auf die Nichtbasisvariablen umgerechneten Zielf unktions- 
koef f izienten £ tab (I) bendtigt wird. 

Existiert zu einer bestimmten Zeile keine Pivotspalte, so ist 

das in der Pivotspaltenindizesma trix durch eine Null (p = 0) 

x z 

angegeben. 

In der Dusenauswahlidentif ikationsnummernmatrix U = Cu ] ist 
fur jede Dusenauswahlidentif ikat ionsnummer I = 1,2, . .., M und 
jede Zeile z = 1,2,..., 6 jeder, in £ abgelegten Koef fizienten- 
matrix die Dusenauswahlidentif ikat ionsnummer derjenigen Dusen- 
auswahl angegeben, die nach dem Basisaus tauschschritt vor- 
liegt, wenn vor dem Basisaustauschschritt die Dusenauswahl I 
vorlag und die Zeile z als Pivotzeile fur den Basisaustausch- 
schritt gewahlt wurde. D.h., wenn zum Zeitpunkt n die Dusen- 
auswahl I[n] vorliegt und ein Basisaustauschschritt mit der 
Pivotzeile z[n] durchgefuhrt wird, kann die Identif ikations- 
nummer der neu entstehenden Dusenauswahl ohne Verwendung von 

A (I) mit I[n + 1] - u _ r - ri direkt aus der Dusenauswahl- 
tab ILnJ z[n] 

identif ikationsnummernmatrix ausgelesen werden. Dies wird 
dadurch moglich, dali durch die Kenntnis des aktuellen Tableaus 
(bzw. seiner Identif ikationsnummer I) und der Wahl der 
Pivotzeile die Pivotspalte implizit festgelegt ist. Zur 
Ermittlung des Folgetableaus ist also die explizite Kenntnis 
der Pivotspalte s[n] unnotig. 

In der Zeilenaustauschindizesmatrix T = ^ij^ ist fur 3 ede 
Zeile z = 1,2, 6 von jeder der in £ abgelegten Koeffizien- 

tenmatrizen der Zeilenindex des Pivotelement s im (sortierten) 
reduzierten Tableau angegeben, das nach dem Basisaustausch- 
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schritt vorliegt, wenn vor dem Basisaustauschschritt die 
Dusenauswahl I vorlag und die Zeile z als Pivotzeile gewahlt 
wurde. 

In der Basisindizesmatrix L = ^^j^ ist fiir j« de Dusenauswahl- 
identif ikationsnummer der Vektor der Indizes der in der Dusen- 
auswahl I enthaltenen Dusen, also der zu I gehorige Basis- 
variablenindizesvektor ( I ) = 1^ angegeben. Die Dusenindizes 
in den Vektoren ^ sind der GroBe nach sortiert. 

Durch den Initialisierungs teil des Verfahrens wird weiterhin 

eine S tart-Dusenaus wahlidentif ikationsnummer I[0] ■ I . und 

start 

die dazugehorige invertierte Basismatrix 

15 — b .start = '-b^start,' ' ^ berei tges tellt . Aus numerischen 

Griinden 

-1 

werden alle Werte in F und B. ^ , die einen sehr kleinen Be- 

— — b. start' 

trag £ , (z.B.: e = 10~ ) haben, zu Null gesetzt. 

20 

Im zweiten Teil des Verfahrens wird zur Realisierung der im * 
Kraf te-Momentenvektor kommandierten Krafte und Momente die 
optimale Diisenans teuerung , d.h. die Ansteuerwerte der 
einzuschaltenden Dusen ermittelt. Nur dieser Teil des 

25 Verfahrens muli den gestellten Echtzeitanf orderungen geniigen. 

Dieser zweite Teil besteht wiederum aus zwei nacheinander 
auszuf uhrenden Teilen. Zuerst wird ein dual zulassiges 
reduziertes Tableau als Starttableau ermittelt. Im zweiten Teil 
wird die erf orderliche Anzahl von Bas i saus tauschschr it ten 

30 durchgef iihrt . 

Ein dual zulassiges reduziertes Starttableau wird wie folgt 

ermittelt: Fiir die Berei tstellung eines dual zulSssigen 

reduzierten Tableaus ist lediglich die Transformation des 

35 aktuellen Kraf te-Momente-Vektors b, auf die fur den Start der 

— kom 

Basisaustauschschritte festgelegte Dusenauswahl I[0] = I 

start 

erforderlich. Dazu dient die schon bereitges tellte invertierte 
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-1 

Basismatrix B b -8tart , und es gilt: 
-1 

-tab. start = B b. start * -konT 

Mit der schon bereitges tellten Dusenauswahlidentif ikat ions- 

nummer fur den Start I [ 0 3 = I t und dem neu berechneten Vek- 

start 

tor ^ tab = —tab start steht ein reduziertes und dual 

zulassiges Starttableau [l[0] b^ w [0]] bereit. 

— tab 

Dann erfolgt die Ausfuhrung der Basisaustauschschritte (zykli- 
scher Verf ahrens teil ) : 

Der folgend beschriebene Verf ahrens teil muB wiederholt solange 
ausgefuhrt werden, bis die optimale Dusenauswahl mit den 
optimalen Ans teuerwerten gefunden wurde Oder festgestellt wird, 
dafi es fiir den aktueli vorliegenden Kraf te-Momenten-Vektor und 
die vorliegende Dusenanordnung keine Losung gibt* 

Da von einem dual zulassigen Simplextableau in Form des redu- 
zierten Starttableaus ausgegangen wird und durch die Basisaus- 
tauschschritte immer wieder nur eine Dusenauswahl erzeugt wer- 
den kann, deren Simplextableau dual zulassig ist, kann bei der 
Durchfuhrung der Basisaustauschschritte immer wieder auf die im 
Initialisierungsteil des Verfahrens bereitgestellten Daten 
insbesondere im Feld der Koef f izientenmatrizen aller optimalen 
und damit gleichzeitig aller dual zulassigen Simplextableaus 
und in der Pivot spal tenindizesma trix zurlickgegri.f f en werden. 

Erster Schritt: 

Finde einen negativen Wert in b K Cn]. n in eckigen Klammern 

t a o 

bezeichnet die Nummer des Iterationsschrittes . Enthalt b^ w [n] 

—tab 

keinen negativen Wert, liegt die Losung des linearen Optimie- 
rungsproblems vor, und es kann der Vektor der Diisenan- 
steuerwerte k Qpt aufgebaut werden. Enthalt £ tab I>] negative 
Werte, dann wird eine dieser Zeilen als Pivotzeile z[n] 
gewahlt. Es ist gunstig, die Zeile des betragsgrdOten negativen 
Elements in b ab Cn] als Pivotzeile zu wahlen. 
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Zweiter Schritt: 

Mit Kenntnis der Diisenauswahlidentif ikat ionsnummer l[n], und 

der im ersten Schritt gewahlten Pivotzeile z[n] wird aus der 

5 Pivotspaltenindizesmatrix P der dazugehorige Pivot spaltenindex 

mit s[n] - p i[ n ] 2 [ n ] ausgelesen. 1st in J? kein Pivotspal- 

tenindex angegeben (p r - r - = 0),mufl eine andere Pivotzeile 

I[n]z[n] 

ausgewahlt -werden. 1st zu keiner der mdglichen Pivotzeilen in 
der Pivotspaltenindizesmatrix eine Pivotspalte indiziert, 
10 existiert keine Losung des linearen Opt imierungsproblems . Sind 

Pivotzeile und Pivotspalte ermittelt, wird zum dritten Schritt 
ubergegangen. 

Im dritten und im vierten Schritt wird der eigentliche Basis- 
15 austauschschritt durch Ermittlung des nachsten reduzierten 

Simplextableaus ausgef iihrt : 
Dritter Schritt: 

Mit Hilfe der Koef f izienten in der Pivotspalte s[n] von A . [n] 

— tab 

wird aus b [n] der neue Basi svariablenwert evektor b^ . 

tab — tab.neu 

20 berechnet. Die bendtigten Elemente des Pivot spal tenvektors a 

aus A b Cn] konnen direkt aus dem Feld der Koef f izienten- 

matrizen aller Optimal tableaus mit a^ . , P -,[n] = f, r -._ P 1 

tab k s[n] k sLnJILnJ 

ausgelesen we r den . 

Die Elemente des neuen Basisvariablenwert evektor s werden wie 
25 folgt berechnet: 

Der neue Wert in der Pivotzeile z[n] wird berechnet nach: 

b = bfttb Z W ^ 

^["l *[«]/[«] 

30 wobei f ^r-i T r n da s entsprechend indizierte Pivotelement 

zLn J sLnJ I[n] 

aus — tab^ n ^ Lst * 

Die neuen Werte in den anderen Zeilen werden berechnet nach: 

3 5 b ubjuu k = b tab k [n]~ b tabMU z[n] • f ks[n] nn] 

wobei : (k = 1(1)6) aber [k * z[n]} 

Da der neuberechnete Basi svariablenwertevektor . zum 

— tab.neu 
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fiktiv angenommenen neuen sortierten Basisvariablen- 

indizesvektor *X fa [n + 1] passen mufl, ist im allgemeinen der 

neuberechnete Basisvariablenwertevektor entsprechend 

umzuordnen. Dazu sind im neuberechneten Basisvariablen- 

wertevektor —tab neu ledi 9 licn die Werte zwischen der aktuellen 

Pivotzeile z[n] und der mit Hilfe von I[n] und z[n] aus der 

Zeilenaustauschindizesmatrix T = [t„. ] auszulesenden Austausch- 

— Iz 

zeile z.^.^utn] « t _ r , r , zyklisch zu tauschen. Die 

tausch ILn] z[n] 2 

Austauschzeile z. . [n] ist die Zeile, wohin der Wert des 

tausch 

neuberechneten Basisvariablenwertevektors plaziert wird, der in 
—tab. neu in der Piv °tzeile zCn] stent. Nach diesem Umordnen des 
neuberechneten Basisvariablenwertevektors liegt der fur den 
folgenden Basisaustauschschritt benotigte Basisvariablen- 
wertevektor b' tab Cn + 1] vor. 

Vierter Schritt: 

Mit Kenntnis von I[n] und der Pivotzeile z[n] wird mit 

I[n + 1] - u r , r die Diisenauswahlident if ikationsnummer 
I L n J zLnJ 

ermittelt, die zum neuberechneten und nach auf steigenden 

Basisvariablenindizes sortierten Basisvariablenwertevek- 

tor £L tab £ n + 1 3 gehort. Der Basisaustauschschritt ist damit » 

vollzogen, das neue reduzierte Tableau Cl[n + 1 ] b^ .[ n + 1]] 

— t a b 

liegt vor, und es kann die Losung des linearen Optimierungs- 
problems mit dem ersten Schritt des zyklischen Verf ahrens teils 
fortgesetzt werden. 

Bricht das Verfahren mit einem Basisvariablenwertevektor, der 

keine negativen Werte mehr enthalt, ab, so liegt der optimale 

6-dimensionale Basisvariablenwertevektor b^_ . = b^ [n ] 

—tab. opt —tab end 

vor, und es kann mit Hilfe von l[n .] aus der Basisindizes- 

end 

matrix L der zum aktuellen Basisvariablenwertevektor gehorende 

Vektor der Basisvariablenindizes mit \b[n ^] = l Tr n 

— end ILn j 

end 

ausgelesen werden. Mit Kenntnis der beiden Vektoren b^ . [n ,] 

— tab end 

und -b^ n end^ w ird der optimale D-dimensionale Ansteuervektor 
fur die Diisen j£ opt aufgebaut. 

Fur die Losung des linearen Optimierungsproblems fur 
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unterschiedliche Kraf te-Momenten-Vektoren sind im allge- 

meinen auch unterschiedlich viele Basisaustauschschrit- 

te n end e tO, 1, n max ^ erf orderlich. Bei Simulationsrech- 

nungen mit einer Diisenanordnung von 12 Diisen wurden in sehr 

seltenen Fallen bis zu n =10 Basisaustauschschrit te durch- 

max 

gefiihrt, obwohl im Durchschnitt nur ca. 3 , 5 Ba s i saus tausch- 
schritte notig waren. Bei geeigneter Wahl des Starttableaus 
[I gtart —tab start^ sind jedoch nur in ca. 0,5 % aller Falle 
mehr als 5 Austauschschritte notwendig. Im allgemeinen bleibt 
es ohne merkliche Konsequenzen fiir die Anwendung, wenn die 
Anzahl der maximal auszuf iihrenden Austauschschritte auf einen 
Wert begrenzt wird, bei dem nahezu alle Falle vollstandig 
bearbeitet werden und in den Fallen, in denen doch noch 
negative Werte im Basisvariablenwertevektor stehen, diese zu 
Null gesetzt werden. Auf diese Weise kann der Worst-Case- 
Rechenauf wand weiter gesenkt werden. Anwendungen dieser 
Verfahrensmoglichkeit bei der Simulation in geschlossenen 
Regelkreisen zeigten keine bemerkbare Verschlechterung der 
Regelgiite oder eine Erhohung des Treibstof f verbrauchs . 

Enthalt der Basisvariablenwertevektor nach der Ausfiihrung der 
maximal erlaubten Anzahl von Bas i sautauschschri tten noch 
negative Werte, ist es prinzipiell auch moglich, dafl keine 
Losung des linearen Optimierungsproblems existiert. Dieser Fall 
kann dann auftreten, wenn eine Diisenanordnung vorliegt, die 
nicht beliebige Krafte und Momente aus dem Bereich aller 
zugelassenenen Kraf te-Momenten-Vektoren erzeugen kann. Damit in 
diesen Fallen keine zu groflen Fehlkrafte und Fehlmomente 
erzeugt werden, kann in Abhangigkeit vom Betrag der negtiven 
Elemente im Basi svariablenwert evektor entschieden werden, ob 
die negativen Elemente zu Null gesetzt werden und der Ansteuer- 
vektor fiir die Diisen aufgebaut wird oder generell keine An- 
steuerung der Diisen erfolgt. 

Kommen bei der Abarbeitung des Simplexalgorithmus entartete 
Optimaltableaus vor, dann ist es theoretisch moglich, daB der 
Algorithmus in einen Endloszyklus lauft und nicht konvergiert. 
In Simulationsrechnungen mit unter schiedlichen Diisenanordnungen 
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ist dieser Fall allerdings noch nie eingetreten und wird 
deshalb auch nicht explizit beriicks ichtigt . Fur den dualen 
Simplexalgorithmus kann man durch die gewahlte Diisenanordnung 
sicherstellen , dafi keine dual entarteten Opt imal tableaus 
exisi:ieren. FUr den priraalen Simplexalgorithmus, also das 
Auftreten primal entarteter Optimaltableaus , ist das so nicht 
moglich. 



Soil bei der Ermittlung der treibs tof f minimalen Triebwerks- 
10 ansteuerung der unterschiedliche spezifische Treibstoff- 

verbrauch , j=1(1)D der einzelnen Triebwerke berucksichtigt 
werden, liegt die Zielfunktion 

j^J k 0Pt J = Min 1 
15 vor. J 

Die Beriicksichtigung der Gewicht skoef f izienten ist in dem 
Verfahren fur den on-line auszuf iihrenden Verf ahrensteil ohne 
Mehraufwand moglich. 

20 

Fallt eine Diise aus, andert sich die Diisenanordnung im Sinne 
dieses Verfahrens. Das Verfahren auf die neue Diisenanordnung 
anzuwenden erfordert, daB alle Daten, die im Initiali- 
sierungsteil des Verfahrens bereitgest ellt wurden, fiir die neue 
25 Diisenanordnung bereitgestellt werden* 

Da sich die neue Diisenanordnung nur durch den Wegfall einer 
Diise von der alten Diisenanordnung unterscheidet , gibt es 
Moglichkeiten, die neuen Daten durch eine Aktualisierung der 

30 alten Daten zu ermitteln. Durch die Bereit s t ellung eines 

erweiterten Feldes von Koef f izientenmatrizen , in dem auch 
diejenigen Koef f izientenmatrizen abgelegt sind, die durch den 
Ausfall von Diisen Koef f izientenmatrizen von Optimaltableaus 
werden, und entsprechend erweiterter AdreBmatrizen, kann der 

35 Aktualisierungsauf wand stark reduziert werden und fiir bestimmte 

Matrizen ganz entfallen. 

SchlieBlich besteht die Moglichkeit, die Daten fiir alle 
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Diisenanordnungen, die durch den Ausfall einer DUse zu erwarten 
sind, vorab bereitzus tellen , und bei Ausfall einer bestimmten 
Diise dann auf die ent sprechenden neuen Daten umzuschalten . Bis 
zum Ausfall einer nachsten Diise miissen fiir die dann zu 
erwartenden Diisenanordnungen die Daten wiederum berei tgestellt 
werden, urn dann durch Umschalten auf den nachsten Dusenausfall 
reagieren zu konnen. 

Auch bei einer Verschiebung des Schwerpunktes des Raumfahr- 
zeuges kommt es zu einer Veranderung der Dtisenanordnung im 
Sinne dieses Verfahrens. Im allgemeinen fiihrt das zu einer 
Veranderung der Systemmat rix. Daher miissen alle von ihr abge- 
leiteten Daten neu bereitgestellt werden. 

Das zuvor beschriebene Verfahren kann zur weiteren Einsparung 
von Rechenauf wand und Speicherplatzbedarf noch modifiziert 
werden : 

- In vorans tehender Beschreibung wurde stets von der Steuerung 
aller 6 Freiheit sgrade des Raumf ahrzeugs ausgegangen. Das 
Verfahren ist auch auf die Steuerung von 1 bis 5 Freiheits- . 
graden anwendbar. Es ist lediglich die Streichung bzw. die 
Nichtberiicksichtigung der den nicht zu steuernden Freiheits- 
graden entsprechenden Zeilen in der Systemmatrix und der von 
ihr abgeleiteten Vektoren und Matrizen erf orderlich . 
Desgleichen ist die Steuerung einer sich standig andernden 
Teilmenge aller 6 Freiheit sgrade moglich. Liegt eine schiefe . 
Diisenanordnung vor, bleiben jedoch die nicht gesteuerten 
Freiheitsgrade im allgemeinen nicht unbeeinf luBt . 

- Miissen am Raumfahrzeug nicht beliebige Krafte und Momente 
ausgeiibt werden, ist der Einsatz einer reduzierten Diisen- 
anordnung moglich. Das vorliegende Verfahren ist auch fiir 
derartige Anordnungen anwendbar. Es sollte nur gewahrlei s tet 
sein, daB die kommandierten Krafte und Momente des Kommandos 
durch die vorliegende Diisenanordnung realisierbar sind. Ist das 
Kommando nicht realisierbar, gibt das Verfahren an, dafl keine 
Losung des linearen Optimierungsproblems existiert. 
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- In Abhangigkeit von der vorliegenden Dusenanordnung kommt es 
vor, daB von den im Peld der Koef f iz lent enmat rizen aller 
Optimaltableaus abgelegten Koef f izientenmatrizen bestimmte 
5 Spalten nie als Pivotspalte ausgewahlt werden. Nach Berechnung 

der Pivotspaltenindizesmatrix P konnen die nicht bendtigten 
Spalten aller Koef f izientenmatrizen einfach ermittelt werden. 
Es ist moglich, ein Feld aller relvanten Pivotspaltenvektoren V 
zu erstellen, in dem nur die Spalten der Koef f izientenmatrizen 

10 aller Optimaltableaus abgelegt sind, die tatsachlich als Pivot- 

spaltenvektoren benotigt werden. Dadurch kann eine signifikante 
Speicherplatzersparnis erreicht werden. Die Adressierung der 
Matrixspalten im Feld aller relevanten Pivotspaltenvektoren 
erfordert keinen erhdhten Rechenauf wand. Jede Spalte hat im 

15 Feld aller relevan-ten Pivotspaltenvektoren ihren Spal tenindex. 

Wird in einer veranderten Pivotspaltenindizesmatrix, der 
Pivot spalteniden-tif ikationsnummernmatrix 0 = ^i^^ nicht der 
Spaltenindex eines Pivot spal tenvektors in einer 
Koef f izientenmatrix angegeben, sondern gleich die 

20 Pivotspaltenidentif ikationsnummer als Spaltenindex r im Feld 

aller relevanten Pivotspaltenvektoren V, dann kann im zweiten 
Schritt des zyklischen Verf ahrens teils rait r = o Ik gleich der 
Spaltenindex im Feld aller relevanten Pivot-spaltenvektoren V = 
[v^^] angegeben werden. Die Ermittlung der Elemente des neuen 

25 Basisvariablenwertevektors im dritten Schritt des zyklischen 

Verfahrensteils andert sich dann wie folgt: 

b = b * ab ffrj ^ 

tabwzln] v z[n]r[n] > 

wobei v z [ n ] r [ n ] das entsprechend indizierte Pivotelement aus 
A tafa [n] ist. Die neuen Werte in den anderen Zeilen werden 
berechnet nach: 

b mbjuu ft = b<ab ft W " b ubMg z[n] ■ V k 
wobei: = 1(1)6) aber (k * z[n\) 
Die weiteren Verf ahrensbestandteile bleiben unverandert. 

- Liegt ein Feld aller relvanten Pivotspaltenvektoren V vor, 
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kann insbesondere bei Dusenanordnungen , die bestiramte 
Symmetrien aufweisen, festgestellt werden, dafi relevante 
Pivotspaltenvektoren unterschiedlicher Optimaltableaus 
numerisch identisch sind. Da das mehrmalige Ablegen gleicher 
5 Pivotspaltenvektoren nicht notwendig ist, kann erneut weiterer 

Speicherplatz eingespart werden. Beim sukzessiven Aufbau des 
Feldes aller relevanten Pivot spaltenvektoren im Initiali- 
sierungsteil des Verfahrens kann iiberpriift werden, ob ein 
numerisch identischer Pivotspaltenvektor zum gerade ermittelten 

10 Pivotspaltenvektor schon in V abgelegt wurde. Ist dies der 

Fall, kann das erneute Ablegen dieses relevanten 
Pivotspaltenvektors entfallen und in der zu V gehorigen 
Pivotspaltenidentif ikationsnummernmat rix () wird die 
Pivotspaltenindentif ikationsnummer ( Spal tenindex) des schon im 

15 Feld aller relevanten Pivotspaltenvektoren abgelegten 

Pivotspaltenvektors eingetragen. 

- Besteht das Erfordernis, den Speicherbedarf des Verfahrens 
noch weiter zu reduzieren, kann eine Verf ahrensvariante 

20 eingesetzt werden, bei der es nicht erforderlich ist, alle 

relevanten Pivotspaltenvektoren, sondern nur noch Reprasen- , 
tantenvektoren im Rechnerspeicher abzulegen. Reprasentanten- 
vektoren sind relevante Pivotspaltenvektoren, aus denen mit nur 
minimalem Rechenauf wand weitere relevante Pivotspaltenvektoren 

25 berechnet werden konnen. 

T 

Jeder Pivotspaltenvektor ji = [a^a^ ### a ^ ### a n ^ gehort zu 
einer bestimmten Basis, z,B. 

X. ( I ) = [X. - „ • • • "X • • • A ] und zu einer 

— b 1 k.11.2 s.i p.n 

30 Nichtbasisvariablen, z.B.X ^ ■ ^ • Wenn a. verschieden von 

nb r x T 

Null ist, laBt sich aus ad^^) = la^a^ ••• a i ••• a^l 



derjenige Pivotspaltenvektor 

a = a(I-,X ) = [a ,a - ••• a ••• a ] mit 

-neu — 2' s neu.1 neu.2 neu.i neu«n 



*5 a = -1 



neud - 57 und a^^k * i) = - a k • a^j 

berechnen, der zur Basis 



20 



WO 98/49058 



PCT/EP98/01595 



Nichtbasisvariablen X _ = X gehort. So kann aus a(I„, ) mit 

no s — i r 

der Angabe des Austauschindexes i der neue Pivotspaltenvektor 
a(I 0 , "k ) berechnet werden. 

— £. S 

Auf diese Weise konnen aus einem z.B. sechsdimensionalen 
Pivotspaltenvektor bis zu sechs weitere Pivot spalt envektoren 
berechnet werden. Fur die Anwendung im Verfahren zur treib- 
stof f minimalen Diis enans teuerung kann man sich jetzt auf das 
Speichern derjenigen relevanten Pivotspal tenvektoren als 
Reprasentantenvektoren beschranken, die erforderlich sind, urn 
alle relevanten Pivotspaltenvektoren wie angegeben zu 
berechnen. 

Fur die Ermittlung einer moglichst kleinen Anzahl von 
Reprasentantenvektoren lassen sich sehr giinstig Genetische 
Algorithmen als Optiomierungsverf ahren einsetzen. 

Fiir eine geeignete Adressierung der Reprasentantenvektoren in 
der Matrix aller Reprasentanten sind statt der 
Pivotspaltenindizes - bzw. 

Pivotspaltenindentif ikationsnummernmatrix einmal eine 
Reprasentantenidentif ikationsnummernmatrix R und weiterhin eine 
Austauschindizesmatrix S, erforderlich. In der 
Reprasentantenidentif ikationsnummernmatrix R wird fiir jede 
Dusenauswahlidentif ikationsnummer und jede Pivotzeile die 
Identif ikationsnummer des Reprasentantenvektor s angegeben, der 
fiir die Berechnung des tatsachlich benotigten 
Pivotspaltenvektors verwendet werden kann. In der 
Austauschindizesmatrix S wird fur jede 

Dusenauswahlidentif ikationsnummer und jede Pivotzeile der 
Zeilenindex des Aus tauschelements angegeben, um aus dem 
Reprasentantenvektor den tatsachlich benotigten Pivotspalten- 
vektor berechnen zu kdnnen. Steht an einer Stelle in der 
Austauschindexmatrix eine Null, bedeutet dies, daB der 
adressierte Reprasentantenvektor unmittelbar der erf orderliche 
Pivotspaltenvektor ist und daher die Umrechnung in den aktuell 
benotigten Pivotspaltenvektor entf a lit . 
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Da stets mit nach Bas isvariablenindizes sortierten 
Pivotspaltenvektoren gearbeitet wird, muB der zunachst aus dem 
Reprasentantenvektor berechnete Pivotspaltenvektor noch 
umsortiert werden (z.B. nach auf steigenden 
5 Basisvariablenindizes) , um benutzt werden zu konnen. Wird die 

Sortierung der Grofle der Basisvariablen vorgenommen, gibt es 
bei einem sechselement igen Vektor genau 30 Moglichkei ten der 
Umsortierung, um die erf or-derliche Ordnung wieder 
herzus.tellen. Bei einer erf orderlichen Umsortierung wird stets 
10 ein Vektorelement an einer Position im Vektor entnommen und auf 

eine andere Position geschrieben. Die Vektorelemente zwischen 
den beiden Positibnen werden ent spre-chend verschoben. 

In einer Sortieridentif ikationsnummernmatrix G = Cg z R ] wird 
15 fur jeden Repra sentant en und jede Aus tauschzeile die Sortier- 

identif ikationsnummer angegeben, nach der der Vektor umgeordnet 
werden muB, um den benotigten Pivotspaltenvektor zu erhalten. 

Aus einem Reprasentantenvektor konnen durch Basi saustausch 
20 und/oder oben beschriebene Umordnung bis zu 7 x 31 - 217 

Vektoren gebildet werden. Werden alle diese Moglichkei ten 
berucksichtigt , kann es aufgrund von Symmetrien in der 
Triebwerksanordnung vorkbmmen, dafl weniger Reprasentanten- 
vektoren erforderlich sind, als wiirden nur die sieben aus einem 
25 Reprasentantenvektor gebildeten relevanten Pivotspaltenvektoren 

in Betracht gezogen, 

Eine weitere Verf ahrensvariante ergibt sich daraus, daB auf den 
Basisaustausch am Reprasentantenvektor verzichtet wird und nur 

30 die 30 Umordnungsmoglichkeiten des Reprasentantenvektors 

berlicksichtigt werden. Auf diese Weise reprasentiert ein 
Reprasentantenvektor 31 Vektoren. Aufgrund von Symmetrien in 
der Triebwerksanordnung konnen aiich auf diese Weise mehrere 
relevante Pivotspaltenvektoren generiert werden. Im Vergleich 

35 zur letzten Verf ahrensvariante werden im allgemeinen jedoch 

mehr Reprasentantenvektoren benotigt, jedoch entfallt 
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gleichzeitig der Rechenauf wand fiir den Basisaustausch am 
Reprasentantenvektor , und es ist die Austauschindizesmatrix £. 
nicht erf orderlich. 
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Patentanspriiche 

1 . Verf ahren zur insgesamt t reibstof f minimalen rechnerge- 

stiitzten Ansteuerung einer beliebig festgelegten Anzahl und 

beliebig festgelegten Anordnung von Diisen an einem 

Raumf ahrzeug, 

dadurch gekennzeichnet , 

- daB ein treibstof f minimaler Anst euervektor (k . ) mit 

opt 

ausschliefllich nicht negativen Ansteuerwerten fiir alle 
vorgegebenen Diisen mit Hilfe des dualen Simplexalgorithmus 
gebildet wird, 

- daB dazu jeder Dusenauswahl aus der Gesamtzahl von Diisen am 

Raumf ahrzeug, die die optimale Dusenauswahl fiir die treib- 
stof fminimale Realisierung eines beliebigen Kommando- 
vektors (b kQm ) darstellt, eine Diisenaus wahliden- 
tif ikationsnummer (I) zugeordnet wird, 

- daB fiir die Ermittlung der optimalen Ansteuerung zunachst 

der Kommandovektor (b, ) der die auszuiibende Kraft 

— kom 

und/oder Momentenkomponenten enthalt, durch Multiplication 

mit einer Transf ormationsmatrix ( ( B j)" 1 ) in einen 

start 

Startvektor (b=b . ) transf ormiert wird, 
— —start 

- daB danach der duale Simplexalgorithmus durchgefuhrt wird, 

bis der Vektor der Ansteuerwerte (b) nur nichtnegative 
Ansteuerwerte enthalt Oder aber die Nichtexistenz einer 
Losung angezeigt wird, 

- und daB dabei der duale Simplexalgorithmus realisiert wird, 

indem vorab alle Daten, die fiir die Abarbeitung des dualen 
Simplexalgorithmus notwendig sind und sich off-line (vor- 
ab) berechnen lassen, ermittelt werden, in Datenfeldern 
im Rechnerspeicher abgelegt werden und beim Einsatz des 
dualen Simplexalgorithmus wahrend des on-line-Bet riebes 
durch geeignete Adressierungen aus dem Rechnerspeicher 
abgerufen und in der Rechnung verwendet werden, wobei nur 
mit einem reduzierten Simplextableau gearbeitet wird, in- 
dem lediglich die DUsenauswahlidentif ikationsnummer und 
der dazugehorige Vektor der Ansteuerwerte (b) gefiihrt 
werden. 
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2. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet , 

daB die Transf ormationsmatrix ((B ^ J" 1 ) die invertierte 

start 

Basismatrix derjenigen Diisen darstellt, die zur Startdusen- 

auswahl gehoren, diese Diisenauswahl fiir einen beliebigen 

Kommandovektor optimal ist, diese Diisenauswahl eine Start- 

Dusenauswahlidentif ikationsnummer (I , ^ ) zugeordnet ist und 

start 3 

der Startvektor (b=b ) die Anst euerwerte dieser Diisen 

start 

enthalt • 



3. Verfahren nach einem der Anspriiche 1 Oder 2, 
dadurch gekennzeichnet, 

- dafl jede Diise am Raumfahrzeug eine Identif ikationsnummer 

erhalt, und die Werte im Vektor der Ans teuerwerte fur eine 
15 bestimmte Diisenauswahl stets in einer festgelegten Ordnung 

anhand der Ident i f ikationsnummern der Dtisen sortiert 
werden , 

- daB bei der Abarbeitung des dualen Simplexalgorithmus nach 

der Ermittlung einer Pivotzeilennummer und mit Kenntnis 
20 der Diisenauswahlidentif ikationsnummer der aktuell vor- 

liegenden Diisenauswahl die entsprechende Pivotspalten- • 
nummer aus einer vorab bereitges tellten Pivotspaltenin- 
dizesmatrix ausgelesen wird, 

- daB das Nichtvorhandensein einer Pivotspaltennummer zu der 
25 vorliegenden Diisenauswahlidentif ikationsnummer und zu 

der gewahlten Pivotzeile durch einen besonderen Eintrag 
in der Pivotspaltenindizesmatrix eindeutig angezeigt wird, 

- und daB die Nichtexi stenz einer Lbsung dadurch angezeigt 

wird, daB zu keiner mdglichen Pivotzeile in der Pivotzei- 
30 lenindizesmatrix eine Pivotspaltennummer angegeben ist. 

4. Verfahren nach einem der Anspriiche 1 bis 3, 
dadurch gekennzeichnet, 

- daB der Pivot spaltenvektor der Koef f izientenmatrix des 

35 aktuellen Simplextableaus , dessen Werte fur die Berechnung 

neuer Ansteuerwerte fur die Diisen einer neuen Diisenauswahl 
benotigt werden, mit Kenntnis der Diisenauswahl vor dem 
Basisaustauschschritt und der gewahlten Pivotspaltennummer 
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aus einem vorab bereitges tellten Feld aller moglichen 
optimalen Koef f izientenmatrizen ausgelesen wird, 

- daB der Vektor der Ansteuerwerte fur die Dusen der DUsen- 

auswahl, die nach dem Basisaustauschschritt vorliegt, 
berechnet wird und die Werte in diesem Vektor entsprechend 
der Zuordnung zu den Dusen und entsprechend der festgeleg- 
ten Ordnung der Identif ika tionsnummern fur die Dusen 
sortiert werden. 

5. Verfahren nach einem der Anspriiche 1 bis 4, 
dadurch gekennzeichnet , 

- daB mit Kenntnis der Dusenauswahlidentif ikationsnummer vor 

dem Basisaustauschschritt und der gewahlten Pivotzeilen- 
nummer die Diisenauswahlident if ikationsnummer der jenigen 
Dtisenauswahl , die nach dem Basisaustauschschritt 
vorliegt , aus einer Dusenauswahlidentif ikationsnummern- 
raatrix ausgelesen wird 

- und dafl nach Vorliegen eines Vektors der Ansteuerwerte, der 

nur noch nicht negative Werte enthalt, mit Kenntnis der 
dazugehdrigen Dusenauswahlidentif ikationsnummer der Vektor 
der Dusenindizes aus einer Basisvariablenindizesmatrix » 
ausgelesen wird und die mit ihrer Dusenidentif ikations- 
nummer indizierten Diisen den ent sprechenden Wert aus dem 
Vektor der Ansteuerwerte (b) erhalten sowie die ubrigen 
Diisen den Ansteuerwert Null zugewiesen bekommen und damit 

der Ansteuervektor (k ) bestimmt ist. 

—opt 

6. Verfahren nach einem der Anspriiche 1 bis 5, 
dadurch gekennzeichnet , 

daB nur diejenigen Spalten von Koef f izientenmatrizen dual 
zulassiger Simplextableaus in einer Matrix aller relevanten 
Pivotspaltenvektoren abgelegt werden, die bei der Abarbeitung 
des Verfahrens tatsa'chlich als Pivotspaltenvektoren ausgewahlt 
werden konnen, 

und daB identische relevante Pivotspaltenvektoren nur einmal in 
dieser Matrix abgelegt werden. 
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7. Verfahren nach Anspruch 6, 
dadurch gekennzeichnet , 

- daB fiir die Adressierung der relevanten Pivotspalten- 

vektoren jedem Pivotspaltenvektor in der Matrix aller 
relevanten Pivotspaltenvektoren eine Pivotspalteniden- 
tif ikationsnummer zugewiesen wird und diese in einer 
Pivot spaltenident if ikationsnummernmatrix der art eingetra- 
gen werden, da0 bei Kenntnis der aktuellen Diisenauswahl- 
identif ikationsnummer und der aktuell gewahlten Pivotzei- 
lennummer die Pivotspaltenidentif ikationsnummer des aktu- 
ell benotigten Pivot spaltenvektors auslesbar ist, 

- dafl das Nichtvorhandensein einer Pivotspalte zu einer 

gewahlten Pivotzeile durch einen besonderen Eintrag in der 
Pivotspaltenidentif ikationsnummernmatrix eindeutig 
angezeigt wird, 

- und dafi die Nichtexistenz einer Losung dadurch bestimmt 

wird, daB zu keiner moglichen Pivotzeile in der Pivot- 
spaltenidentif ikationsnummernmatrix eine Pivots pal ten- 
identif ikationsnummer angegeben ist. 

8. Verfahren nach einem der Anspriiche 3 bis 7, 
dadurch gekennzeichnet , 

dafi statt der relevanten Pivot spaltenvektoren nur bestimmte 
Pivotspaltenvektoren als Reprasentantenvektoren im 
Rechnerspeicher abgelegt werden, wobei die 
Reprasentantenvektoren eine mdglichst kleine Anzahl von 
Pivotspaltenvektoren darstellen, aus denen die tatsachlich 
benotigten relevanten Pivotspaltenvektoren ermittelt werden, 
und dafi jeder Reprasentantenvektor eine 

Reprasentantenident if ikationsnummer zugewiesen bekommt und zur 
Adressierung der Reprasentantenvektoren eine 

Reprasentantenidentif ikationsnummernmatrix verwendet wird, in 
der fiir jede Dusenauswahlident if ikationsnummer und jede 
Pivotzeilennummer die Reprasentantenident if ikationsnummer 
desjenigen Reprasentantenvektors eingetragen ist, der fiir die 
Berechnung des erf orderlichen Pivot spalt envektors verwendet 
wird. 
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9. Verfahren nach Anspruch 8, 
dadurch gekennzeichnet, 

daB in einer Austauschmatrix fur jede Dusenauswahl- 
identif ikationsnummer und jede Pivotzeilennummer der 
Zeilenindex der auszutauschenden Basisvariablen als 
Austauschelement angegeben ist, urn aus dem Reprasentantenvektor 
den benotigten Pivot spaltenvektor zu berechnen. 

10. Verfahren nach einem der Anspriiche 1 bis 9, 
dadurch gekennzeichnet , 

daB flir die Moglichkeit der Veranderung der Diisenanordnung 
durch zu erwartende Diisenausf alle die entsprechenden Daten 
vorab berechnet und gespeichert werden und nach Feststellen 
eines Dusenausf alls automatisch durch den Rechner auf Verwen- 
dung der entsprechenden neuen Daten umgeschaltet wird. 

11. Verfahren nach Anspruch 10, 
dadurch gekennzeichnet, 

daB die nach einem Dusenausfall erf orderlichen neuen Daten 
durch eine ent sprechende Korrektur und Aktuali s ierung der bis 
zum Dusenausfall verwendeten Daten bereitgestellt werden. 

12. Verfahren nach einem der Anspriiche 1 bis 11, 
dadurch gekennzeichnet, 

daB das Verfahren unabhangig von der Dimension des Kommando- 
vektors (*i kom ) angewendet wird und dabei fur die Berechnung der 
vorab bereitgestellten Daten die nichtauszuiibenden Komponenten 
aus der vollstandigen Systemmatrix der Diisenanordnung 
eliminiert werden. 
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